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in continuous variable quantum systems. On one hand we derive a Bell-like 
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reconstructed state, thus giving the possibility to use standard Bell's inequalities. 
Illuminating examples are also discussed. 
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1. Introduction 

In their famous paper of 1935 0, Einstein, Podolsky, and Rosen (EPR) questioned the 
completeness of quantum mechanics (QM). Their argument was based on the premise 
of no action at a distance (locahty) and reaUsm. Several years later Bell showed 
that the prediction of QM are incompatible with the premises of local realism (or local 
hidden variable theories). 

Experiments J3 based on Bell's result support QM, indicating the failure of local 
hidden variable theories. Quantitative tests used mainly Bohm's (dichotomic) version 
U of the EPR entangled states instead of the original EPR states with continuous 
degrees of freedom. 

In recent years, systems with continuous variable has attracted much atttention in 
connection with the burgeoning field of quantum information pj. In such a field EPR 
entanglement and quantum nonlocality are of practical importance. Nevertheless, the 
generalization of Bell's inequalities to quantum systems with continuous variabbles 
still remains a challenging issue. 

In continuous variable bipartite quantum systems, EPR aspects can arise when 
trying to infer quadratures of one subsystem from those of the other . These can 
be tested by exploiting the Heisenberg uncertainty principle Such approach has 
been employed in Refs.[3IH]. Few other theoretical proposals use the field quadrature 
measurements |9| 1 1 01 111) , while a phase space approach to quantum nonlocality has 
been developed ^21 and tested [Oj . 

On the other hand quantum tomography |14| provides a useful tool to reach a 
complete state knowledge. Such knowledge could then be used for nonlocal tests. 
Some aspects of EPR problem have been considered in tomographic approach in ■ 

Here, we propose to use the tomographic data to do nonlocal tests. At first 
instance, we shall develop an inequality based on the tomograms. However, this will 
result quite difficult to violate. A state does not have to violate all possible Bell's 
inequalities to be considered quantum nonlocal; rather, a given state is nonlocal when 
it violates any Bell's inequality Thus the degree of quantum nonlocality that 

we can uncover crucially depends not only on the given quantum state, but also on 
the Bell operator 17 . Then, we shall introduce pseudospin operators whose statistics 
can be inferred from the data characterizing the reconstructed state, thus giving the 
possibility to use standard Bell's inequalities. We shall also present some illuminating 
examples. 

2. Quantum tomography revisited 

A quantum state is described by the density operator p or by any its phase space 
representation like the Wigner function W{q,p), where q, p are canonical conjugate 
variables. Quantum tomography is a technique which allows to recover the quantum 
state from a set of measured probability distributions. The latter, sometimes called 
tomograms, are line projections of the Wigner function. Such projections can be 
parametrized through a symplectic transform |18| . Practically, the tomograms can be 
expressed as 



where fj,, v € Tl, are the parameters of the symplectic transformation and X is a 
stochastic variable. It represents the random outcomes of the measurement of the 




(1) 
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observable X — fj,q + i/p, and w results the probability distribution associated to such 
observable. It is worth noting that the above symplectic transformation, once thought 
as composition of rotations and squeezing, can be realized in optical systems I19j . A 
remarkable property of the tomograms ^ is their homogeneity. 

To get the phase space picture of the quantum state one has to invert the relation 

obtaining 

Wiq,p)^ j'w{X,^i,i^)e-""^~"'P+'^ dXd^idi^. (2) 

This simply results an inverse Fourier transform. Neverthless, it becomes a more 
involved inverse transform, namely an inverse Radon transform, when reducing the 
symplectic transformation to a mere rotation, i.e., = cos 9, v = sin 9. This is the 
case of optical homodyne tomography (OHT). 

By expanding the density operator in terms of a complete set of operators, it is 
also be possible to directly relate the quantum state to the tomograms, that is 

p — ^ w{X, ii,h')K{X, ji.v) dX dpidv ^ (3) 

where the kernel operator takes the form 

K{X,fi,iy) = -^X^ cxpi—iXX + i\'^^\exp{iXp.q) exp{iXiyp) . (4) 
zvr \ 2 / 

Here, A can be set equal to unity; this freedom reflects the overcompleteness of the 

information obtainable by means of all possible tomograms 

In OHT it is preferable to use Eq. Q instead of Eq. Q since the kernel in 

the latter case is unbounded, thus requiring the use of an artificial cutoff to directly 

sample the Wigner function. On the contrary, Eq. (j2Jl allows a direct sample of the 

density matrix elements on the basis where K results bounded j2(J| . 

3. Bell-like inequality for tomograms 

The above arguments can be generalized to the case of a bipartite two-mode system 
|21| . In such a case the tomograms will be given by 

w{Xi,ni,iyi ,X2,fi2,i^2) = J dqidpi J W{qi,pi,q2,P2) ^ 

X S{Xi - gi^i -piVi)S{X2 - (72^*2 -P2V2)dq2dp2 , (5) 

where /xi, z^i, /i2, V2 are parameters. Let us restrict to the case of OHT, and suppose 
to perform homodyne detection in both subsystems 1 and 2. Then, p,i = cos 6*1, 
Vi = sin^i, pL2 = cos 02, 1^2 = sin 02, where 0i, 02 are the rotation angles related 
to the local oscillators phase. As a consequence we have /zi, i^i, X2, ^2, ^^2) 

01,^2, 02). 

We now classify the results of the measurements to be + if the quadrature result 
X is greater than zero, and — otherwise Then, we can construct the following 
probabilities from the tomograms 

CO oc 

u;++(0i,02)= J dXi Jw{Xi,9i,X2,92)dX2, (6) 









w+^i9i,92)= J dXi J w{Xi,9i,X2,92)dX2, (7) 
-00 
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oo 

W-+(fc'l 



,02)= J dXi Jw{Xi,9i,X2,e2)dX2, (8) 

-00 


,62)= J dXi j w{Xi,ei,X2,B2)dX2. (9) 



—00 —00 



Then, the Bell's inequalities can be written in terms of the above probabilities. In 
particular, the CHSH inequality (551 can be written as 

B = \E{ei,e2) + E{e^,e'^) + E{e[,e2) - E{e[,B'^)\ <2 , (lo) 

where now 

£^(^1,^2) = w++{ei,e2)-w+.{ei,92)-w.+{ei,92)+w..{eue2)-{ii) 



4. Reconstructed data and pseudo-Spin operators 

As discussed in Section 2, measuring tomograms allows to statistically sample the 
density matrix elements at least in some basis, e.g. Fock basis. This represents a 
complete knowledge of the state, also for a bipartite system. Then, such information, 
say the density matrix elements in Fock basis, can be used to derive the statistics of 
any hermitian operator, even if this is not directly measurable. The price one has to 
pay is the large amount of data needs to collect. 

Let us now see how these arguments can be fruitful applied to the problem of 
quantum nonlocality. We introduce the local pseudo-spin operators |23j 
00 

Si''^ = ^ |2n>,(2n + l| + |2n+l>,(2n| , 

n=0 

Sv'^= -*f;(|2n),(2n + l|-|27i + l),(2n|), (12) 

00 

;^(-)"N,(ni, 

where j = 1,2 and \n)j are Fock states. The operators H12|l obey the commutation 
relation of spin-^ algebra, namely 

= 2 i 5jj' Exyz S^z'^ , (13) 

with Exyz the totally antisymmetric tensor. Spin tomography approach was considered 
e.g. in '5T. 

In terms of the operators H12[) the CHSH inequality |52] reads in its standard 
form, that is 

S= |£:(u,v)+£:(u',v)+f(u,v')-f(u',v')| <2, (14) 

with 

£(u,v)=((u.sW)(v.S(2))). (15) 

Here u, u', v, v' are unit vectors in R^, while S'-') = {Sx \ Sy \ si^'') and the dot 
indicates the ordinary scalar product in R^. Furthermore, the angle brackets in Eq. 



5. 



0) q(/) 
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H15|l denote the average over the density matrix elements assumed available from 
tomographic reconstruction. 

Let us now study the possible violations of inequality (|14|l for the two preceeding 
examples. For simplicity, we consider the directions u, u', v, v' coplanar, in the plane 
X — z, so that 

£(u, v) ^ ^ [sP cos 0u + Si^'' sin Ou) (sf ^ cos 9., + S^^'' sin 0^) ^ , (16) 
where 6'„ {9^) is the relative angle betweeen u (v) and z. 



5. Examples 



Let us going to study the possible violations of the inequalities (|10|) and H15|) in three 
paradigmatic cases. 

5.1. Example A 

We first consider the two-mode squeezed vacuum state which can be generated in the 
non-degenerate optical parametric amplifier |25) 



1^) = Vl-A2£A'»i|n)2 



(17) 



n=0 



where A = tanh s and s is the squeezing parameter. Such state can also be described 
by the following Wigner function 

4 



W{qi,pi,q2,P2) = -^expj-e [(<7i - 92)^ + (pi + ^2)^] 
- e'^[{q^+q2f + {pi-p2r]). 



Notice that for s ^ 00 the state (I18|l becomes precisely the EPR state. 
Inserting Eq. (|18(l into Eq. © we obtain 

2 

TT 



where 



wiXi,9i,X2,92) ^ -Nexp(-2aX?; - laXl - 40X1X3) , 
cosh(2s) 

~ cosh2(2s)-sinh^(2s)cos2(6li-h6l2) ' 
^ _ sinh(2s) cos(6'i + 6*2) 

~ cosh2(2s)-sinh^(2s)cos2(6li+6l2) ' 
X = Va^ - 52 . 

In this case, the tomograms depend only on the sum of the parameters. 
Using the position of Eqs. © we obtain 



(18) 



(19) 



(20) 



u;++(0i,02) =«^--(^i,e2) = 
w+-(9x,92)^w-+{9^,92) = 



1 

47r 
1 

47r 



TT — 2 arctan | — 



TT + 2 arctan 



(21) 



(22) 



where the signs in front of tt refers to 7r/2 < (^i + 92) < 7r/2, and should be adjusted 
accordingly to the range of {9i +92). 
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(a) The quantity ui-)_+(9i, 62) of Eq. 

is shown vs {9i + 62) for dif- 
ferent values of A. The quantity 
{61,62) has identical behav- 
ior with an offset of tt in {61+62). 



(b) The quantity B of Eq. is 
shown vs 6u for different values 
of A. The other parameters are 
6/„ = 7r/4, 6»„/ = -it/2 and 
6/„, = -7r/4. 



Figure 1. Example A (dashed line A = 0.20, dotted line A = 0.54, solid line 
A = 0.96) 



In Fig. 1(a) we show the behavior of the functions H21(l - (|22|l vs 9i + 92 for different 
values of A. As we can see they never exceed the value 1/2, thus the quantities E, 
in Eq. Ullfl . will be always confined between 1 and —1. In turns, the inequality (|l()|l 
will be never violated. This results seems in agreement with Bell's conclusions stating 
that the original EPR state does not exhibit nonlocality because its Wigner function 
is positive everywhere, and as such allows a hidden variable description ^ . However, 
we shall see that is simply a matter of the type of measurements. 

For the state l(T7|) . Eq. lfTC|l becomes 

2A 

£(u,v) = cos cos 6*1, H — sin 9^ sin 9y . (23) 

1 -f A^ 



Then, in Fig. [T^b) we show B for different values of A. In this case the EPR state show 
its nonlocality. Indeed a maximal violation of inequality H14|) is achieved for A ^ 1 
(s oo). 

5.2. Example B 

As second example we consider the entangled superposition 

|^) = -^(|0)i|0)2 + |n)i|n)2) . (24) 
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Since |0) and |n) are orthogonal for n 0, the state H24|) resembles a spin triplet state. 
The corresponding Wigner function is 

1/2" „ „ 

W{qi,pi,q2,P2) = 2:^(1 + -^{qi - ipiT{q2 - ip2f+ 

+ ^{qi+ipiTiq2 + iP2T + (25) 



L„i2ql + 2pl)L„{2qi + 2pl)\e 



2222 

-91-92 -Pi -P2 



where L„ indicates the Laguerre polynomial of order n. By means of Eqs. ^ and 
H25|l we get the tomogram of the state under consideration 

1 / HJXi)HJX2) 
wiXue^,X2,e2) = 2^'i^, '^ osn(gi-02)+ 



+ ^^-^2' y^t-^i. (26) 



HliX^)H^jX2) Y_x?-x? 
22»(n!) 

where iJ„ indicates the Hermite polynomial of order n. Using the position of Eqs. © 
we obtain probabilities w±±{9i,92) 

1 (0) 

w++iei,92)=w--{0i,e2) ^ - + "'J , cosn(^i+g2), (27) 

1 (0) 

w;+_(0i,02) = w^+{9i,d2) = - - "'^ , cosn(gi+g2). (28) 

4 7r2"n! 



In Fig. 2(a) we show the behavior of the w++ vs 61 + 6*2 for different values of 
n. Also in this case they never exceed the value 1/2, thus the quantities E, in Eq. 

will be always confined between 1 and —1, and the Bell's inequality (|10|l will be 
satisfied. 

For the state (f^ . Eq. ifT^ becomes 

£(u,v)^{ '^"^(f"-^;) " = (29) 
^ ^ \ cos6'„ cos 6*1, n > 1. ^ ^ 

Then, in Fig. |2(b)| we show B for the case of n = 1. The nonclassical character of the 
entangled state 1)24}! becomes now manifest. 

5.3. Example C 

As a third, and last, example we consider a pair-coherent state namely 

2-K 

|?A)=iV j \re''P)\re-''^) dif, (30) 


where norm N is 

2 

^ = (31) 

27rVio(2^ 

with /„ the modified Besssel function of order n. The state (13011 represents two 
coherent states, of the same amplitude r G R, having a well defined phase relation 
although their own phase is random. 
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(a) The quantity ui-|_+{9i, 62) of Eq. 
(771 is shown vs (6*1 + 82) for 
different values of n (dashed hne 
n = 1, dotted line n = 3, solid 
line n = 5). 



(b) The quantity B derived from 
Eq. I29i is shown vs 9u for case 
of n = 1. The other parameters 
are 6y = 0, 6„/ = n and 6^1 = 
n/2. 



Figure 2. Example B 



The Wigner function corresponding to the state (|30|l is 

W{qi,pi,q2,P2) = expi-qf - pj - q2 - pI) X 

2-K 2-K 

X J Jexp(V2r(qi-ipi)(e''^ + e-''^') + y2r((Z2-ip2)(e-''^ + 


X exp^r^[2cos((^ — ip') — cos(2i^) — cos{2^p' y\^d^p dip' . 
The tomogram of the state (|3U|I is 

01,^2, 02) = — e-^- 01,^2, 02)1 e-^l-^^ 

TT 

where the integral /(Xi, 0i, X2, 62) reads 

/(Xi, 01, X2, 02) = Jexp [ - ^ (e2»(^-«i) + e-2'(^+^=)) + 
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Using the procedure shown in the Appendix together with some algebra, we can obtain 
from Eq. (|33|l the foUowing expressions for the probabihties w±±{9i,62): 

w±±{ei,92) ^ ^e'^"^' J J cxp(2r^ cos(v3i + (^2)) x 


X ^1 ip erf ^-^ ^e'^vi-vo) ^_ g»(v2+¥'o) j jj y. (35) 
X =F erf J-^ ^^-tivi+vo) ^ g-»(¥^2-¥'o)jjj , 

where erf denotes the error function and ipQ — (6*1 + 92)/ 2. Notice that the first 
(second) ± on the Lh.s. corresponds to the first (second) =F under the integral on the 
r.h.s. 

The probabilities (|35|l are no longer bounded between the values {0, 1/2}, thus, 
they can lead to violation of Bell's inequalities. In particular, we show in Fig. 3(a) the 
dependence of B from r revealing violation of the CHSH inequality in a small range 
of the amplitude r. This result recall that of Ref. i9| where the state l|5n|l was shown 
to violate another Bell's inequality. 

For the state jSUI) . Eq. ifTCjl becomes 

(Jo(27'^) \ 
1 - T sing^ sin6',, , (36) 

where J„ denotes the Bessel function of order n. Then, in Fig. 3(b) we show B for 
r — 1.05. The violation of CHSH inequality becomes now more pronounced. 



6. Conclusions 



In conclusions we have presented a tomographic approach to the quantum nonlocality 
of a bipartite continuous variable system. At first instance we have proposed a rough 
use of tomograms for testing violation of local hidden variable theory. In such a case 
one does not really need of a complete set of tomograms. However, the presence of 
purely quantum correlations are not easily recognizable. As matter of fact, the EPR 
state does not show the nonlocal character, while a pair-coherent state does. 

Then, we have proposed to exploit the tomographically reconstructed state to get 
the statistics of pseudo-spin operators. This, practically, maps a continuous variable 
system into a discrete variable system making possible the use of standard Bell's 
inequalities. This seeems a more powerful approach, but the price one ought to pay 
is the large amount of data (a complete set of tomograms) needed to reconstruct the 
state. In such a case the EPR state completely shows its nonlocal character. Also the 
pair-coherent state evidenciates much larger violations |27j . 

Our result opens possibilities for testing QM against local hidden variable theories 
using very efficient detection methods like homodyne detection pS], 

Finally, the present approach could be straightforward extended to multipartite 
systems, e.g. making possible test of Grenbergcr-Hornc-Zeilinger theorem |2H| in 
continuous variables version j29| . Moreover, it could be useful to characterize quantum 
states used in information processing [5]. 
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(a) The quantity B derived from 
Eq. 1351 is shown vs r for 
parameters 9u = vr/2, By = 
-7r/4, e^, = 0, e,, - 37r/4. 



(b) The quantity B derived from 
Eq. is shown vs du for 

parameters r = 1.05, = 0, 
= TT, = 7r/2. 



Figure 3. Example C 
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Appendix 

By refering to Eq. (|34|l we note that the exponent of the imaginary argument is 
periodic, and since the integral is taken over a period of the exponent we can shift the 
angle argument as (/? — (0i — 6'2)/2 ^ so the integral takes the form 

2 

I{Xi,9i,X2, 62) = Jexp [ - L.e-'^'V'^ (e^^'^ + e'^^'^) + 



where we have set 



Then, the integral can be rewritten as 

2-n 



I{Xi,9i,X2,92) - J exp 



(37) 
(38) 
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X exp --I^ — - + \/2ae-"^X2 dip, (39) 



where we have set 

a = re-'f° . (40) 

Expanding each exponent in the above integral into a series of Hermite polynomials, 
and taking into account the equality 

27r 

Je*("-'")'^d^^27r^„,™, (41) 



we obtain the following expression 

^0 2»(n!)2 
Let us now consider two generic series 



oo 



/(Xi) = Yj /"^l 5(3^2) = 2 grnX^ , (43) 

n— m— 

whose product yields 

00 

h{xi,X2) = fixi)g{x2) ^ ^ fngrnX"x^ . (44) 

n,m— 

How to express the sum of diagonal elements 

oo 

s{xi,X2) = fngnX'lx^ (45) 

in terms of the h{xi,X2) ? The answer is found by taking into account the equality 
H41|l . and it reads 

271- 

s{xi,X2) = ^ J h{xi^'^,X2e-''^) d^. (46) 



The sum in right-hand side of Eq. (|42|l is the sum of diagonal terms of the product of 
two following series 

/,(«) = g^^(^f fc.1,2. (47) 

If we apply the formula (|46() to this sum we obtain the representation H37() for 
it which we have started from. To obtain probabilities w++ etc. we represent 
\I(Xi,ei,X2,e2)\ in the following way: 

27r 27r 

|/|'- J J hiae^nfiic.e^nf2iae'''nf2ic^e^ndip^dip2. (48) 



To calculate probability w++ we need to take the following integral: 

27r 27r 

j j |/(Xi,0i,X2,e2)|'e-^i-^^'d^id(^2, (49) 



and analogously for other probabilities. Swap integration order in these integrals we 
obtain the expression 
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